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We study the fundamental limit on single-photon indistinguishability imposed by decoherence due 
to phonon interactions in semiconductor quantum dot-cavity QED systems. Employing an exact 
diagonalization approach we find large differences compared to standard methods. An important 
finding is that short-time non-Markovian effects limit the maximal attainable indistinguishability. 
The results are explained using a polariton picture that yields valuable insight into the phonon- 
induced dephasing dynamics. 
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The study of the coherence properties of single pho- 
tons emitted from semiconductor cavity QED (cQED) 
systems is important for applications in quantum infor- 
mation technology JL ^^nd provides insight into the funda- 
mental decoherence effects induced by the environment. 
For all-solid-state cQED systems, such as a quantum dot 
(QD) embedded in a photonic crystal cavity [5] [Fig. 
RTa)] or a micropillar cavity [3], the main decoherence 
mechanism at low temperatures is the electron-phonon 
interaction [IHB] , as many recent studies show [21 l7l-fT4] . 

Decoherence limits the degree of indistinguishability of 
single photons emitted from cQED systems [Fig. fl|a)], 
thus diminishing their applicability for scalable linear 
optical quantum computing [T], where an all-solid-state 
single-photon source is a key element. Furthermore, re- 
cent experimental results [H |H1 EH EH [H] necessitate a 
departure from the well understood paradigms of atomic 
cQED, since the strong interaction with reservoirs in the 
solid state calls for new basic models and physical inter- 
pretations. A better understanding of phonon-induced 
decoherence thus leads to insight into the fundamental 
physics of nanostructured solids, and can help ushering 
novel quantum technological devices. However, thus far 
only little attention has been given to the influence of 
phonon interactions on the indistinguishability. Only few 
experiments have been reported [31 [TSl El] and previ- 
ous theoretical studies have employed a Markovian pure 
dephasing approximation [T7l - f24] or phenomenological 
descriptions of finite- memory dephasing processes [25] . 
none of them treating the phonon interaction microscop- 
ically while taking into account the cavity. 

In this Letter we show that the non-Markovian na- 
ture of the phonon reservoir has a large effect on single- 
photon indistinguishability: short-time virtual processes 
occurring on time scales much shorter than a typical 
"dephasing time", must be considered. Also, it is es- 
sential to treat the phonon interaction microscopically 



and on equal footing with the electron-photon interac- 
tion. The analysis is based on an exact diagonalization 
(ED) technique, retaining the inherent non-Markovian 
nature of the phonon interaction to all orders in the 
phonon coupling. Our findings are contrasted to stan- 
dard approximate approaches for including phonon inter- 
actions pLSj 26-29 , namely second order expansions and 
phenomenological pure dephasing descriptions. Figures 
lljc) and (d) show such a comparison. The deviations 
between the approximate and the ED results are signifi- 
cant, demonstrating that memory and back-action effects 
in the reservoir cannot be neglected, as in Markovian ap- 
proaches. 
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FIG. 1. (a) Schematic of a coupled QD-cavity system inter- 
acting with longitudinal acoustical phonons and a photonic 
crystal cavity with an embedded QD. (b) Two-time plane 
with the time diagonal, r = 0, and phonon reservoir corre- 
lation time Tcorr- The shaded region shows the extent of the 
short-time regime, (c), (d) Calculated indistinguishability as 
a function of the QD-cavity coupling strength for light emit- 
ted from the QD and the cavity. Parameters: F = 0.5 ns~^, 
hn = 125 ^eV, ftA = 27.78 neV, and hj = 0.85 peV. 



To calculate the indistinguishability, we model the cel- 
ebrated Hong-Ou-Mandel experiment [15l [30] where two- 
time correlation functions for the photon operator a need 
to be considered [IH], {a){t')a{t)) = {a!^ {t+T)a{t)) . These 
are defined in the plane spanned by t and t' [Fig. [l|b)], 
whereas one-time correlation functions, {a) {t)a{t)) , re- 
side on the time diagonal (r = 0). Consequently, one- 
time functions only experience short-time non-Markovian 
phonon effects within the phonon reservoir correlation 
time, Tcorr [Fig. [l[b)], after the initial excitation. In con- 
trast, the two-time function implies that a photon is re- 
moved, a{t), at each instant t and added again, a\t + T)^ 
an instant t later. This results in a "continuous excita- 
tion" of short-time transients, illustrated as the band sur- 
rounding the time-diagonal in Fig. Iljb). Short-time non- 
Markovian effects thus play an important role through- 
out the entire lifetime of the excitation, strongly affecting 
physical quantities derived from two-time functions, such 
as the indistinguishability which we will demonstrate. 

Theory. - To model a QD coupled to a cavity in- 
teracting with longitudinal acoustical phonons, we em- 
ploy the Jaynes-Cummings model including the electron- 
phonon interaction [3 H [101 [H] [Fig. Qa)]- We fol- 
low Hohenester [32] and employ a set of effective phonon 
modes. Expanding the QD-cavity system in the basis 
{|1) = |e,n = 0), |2) = |g, n = 1), |3) = |g, n = 0)}, where 
n is the cavity photon number, the total Hamiltonian be- 
comes 



H = H 
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where CTy = |j)(j|. The second term describes the phonon 
interaction with the QD, p denotes the effective phonon 
modes with bosonic operators bp and &|,, and Mp is the 
effective phonon matrix element 03] ■ The last term is the 
free phonon Hamiltonian with LOp denoting the frequency 
of mode p [31] . The Jaynes-Cummings model is 



ffjc = ^Ao-ii + hg{(Ti2 + o'2i), 



(2) 



where A = u^g — Wcav is the QD-cavity detuning with uj^g 
and Wcav being the QD and cavity transition frequencies, 
respectively, and g is the QD-cavity coupling strength. 

The system dynamics is obtained by employing the re- 
duced density matrix formalism including Lindblad de- 
cay terms [351 [3S| • The Master equation for the density 
matrix is 



dtp{t) = -in-^[H,p{t)] + Si^p{t), 



(3) 



where the Lindblad terms are Shp{t) = {L^ {0-32} + 
ir {o'3i})p(i). The rate n describes the escape of cav- 
ity photons, related to the Q-factor as ac = Q/wcav, and 
the rate F describes the decay of the QD in the absence 
of the cavity. The Lindblad operator is L^{0}p(t) = 
-2 [O^Op{t) + p{t)0'^0 - 20/9(t)Ot] . Importantly, the 



electron-phonon interaction is here included in the uni- 
tary part of the Master equation and not via approximate 
scattering terms. Thus, the electron-phonon interaction 
is treated on equal footing with the electron-photon inter- 
action, ensuring a rigorous inclusion of all non-Markovian 
phonon effects. To solve Eq. ([3]), we expand the phonon 
operators in a multi-phonon Fock state basis |37j . and 
propagate the equations numerically, providing an exact 
diagonalization of the coupled QD-cavity-phonon system. 
In the limit of 5 — >■ 0, our model becomes the exactly solv- 
able independent boson model (IBM) [35tiiT] . The IBM 
can thus extrapolate the ED results to this limit, where 
the ED becomes cumbersome due to long numerical in- 
tegration times. 

We compare our simulations with two standard ap- 
proaches. The first treats the electron-phonon inter- 
action to second order using the time-convolutionless 
method (TCL) [36l SlIlS]. Here the density operator 
is pit) = Trphon [p(i)], where the phonons are traced out 
and hence treated as a thermal reservoir. Within this 
approximation, the Master equation becomes 

dt-p{t) - ^ih-^ [Hjc, P{t)] + STCL{t)pit) + 5Lp(i), (4) 

where the effects of the phonons are contained in S'tcl (t) 
[TUl H3] . The second approach is a Markovian Lindblad 
description of the pure dephasing processes, equivalent to 
the TCL for a memory-less phonon reservoir. The conse- 
quence is the replacement of S'tcl (i) with the Lindblad 
operator L2.y {cru}, where 7 is the pure dephasing rate 
[25y^ . The rate 7 will be chosen to provide a reasonable 
fit to the ED. 

To calculate two-time functions we invoke the Quan- 
tum Regression Theorem (QRT) [35], which does not 
imply any approximations in describing non-Markovian 
phonon effects in the ED approach. Applying the QRT 
to the TCL density matrix in Eq. Q, however, requires 
more subtle considerations. The TCL results in time- 
dependent scattering rates [TOl [42] and thus the substi- 
tution 7tcl(0 ~^ 7tcl(t) might be expected to include 
non-Markovian effects in the QRT. This is however not 
the case, as recently shown by Goan et al. [101 [JT] . How- 
ever, in the long-time limit, S^chit — ^ 00), the QRT 
becomes a consistent approximation. 

The indistinguishability of the emitted photons is 
quantified as the normalized number of coincidence 
events at the HOM output detectors and is calculated 
as [IE1121 



/ = 



J„^dtJ^°^dr\{AHt + r)Am^ 
/r di/o°° dr{AHt + T)Ait + r))(At(i)A(t)) ' 



(5) 



where A is either the photon operator a = \n = 0){n = 
1| — 1732 for light emitted from the cavity, or the QD op- 
erator a_ — \g){e\ = 0-3 1 for light emitted from the QD. 
We assume a 1-photon basis and second-order contribu- 
tions are therefore absent. 



For all simulations, except when varying the detuning, 
a QD-cavity detuning equal to the polaron shift [lOj has 
been used, ?iA = /lApoi ~ 27.78 ^eV, corresponding to 
an effective QD-cavity detuning close to zero. To investi- 
gate the fundamental limits on indistinguishability set by 
the vacuum phonon bath we set T = K. We emphasize 
that neither our model nor approach are limited to zero 
temperature, however, finite temperatures significantly 
increase the computational effort. The initial condition 
is the QD in the excited state, with both the photon and 
phonon fields in their ground states, corresponding to the 
experimental situation of excitation of the system with a 
short optical pulse, usually employed in measurements of 
the indistinguishability. We neglect effects such as timing 
jitter and nearby fluctuating charges, as these depend on 
the excitation mechanism and can be avoided. 

Dependence of QD-cavity coupling strength. - Figures 
W[c) and (d) show the indistinguishability for light emit- 
ted from the cavity and QD as a function of the QD- 
cavity coupling strength g. The ED and IBM results [44] 
differ quantitatively and qualitatively from both the TCL 
and Lindblad results. For small but increasing g the in- 
distinguishability remains constant in the ED and IBM 
approaches for the QD, which is also expected for the 
cavity, where the IBM does not apply. The TCL may 
predict an indistinguishability above unity [Fig. Iljc)], 
which is unphysical and a well-known issue associated 
with this method 03]. Further increasing the QD-cavity 
coupling, g, the indistinguishability decreases for the ED 
and TCL, whereas the Lindblad theory predicts satu- 
ration. The surprisingly large deviations between the 
ED and common approaches are important, especially 
for applications with strict requirements on the indistin- 
guishability, e.g., in implementations of linear quantum 
computing protocols [Tj . 

The Lindblad theory can qualitatively be interpreted 
using [m [22] / = reff/(reff + 27), where Feff is the 
effective decay rate of the QD. For large detuning or 
small QD-cavity coupling, compared to the loss rates, 
one obtains EH]: Feff = F + 2g^-ftot/{A^ + 7?ot), 7tot = 
1/2(k -I- F -|- 27). The expression predicts an initial in- 
crease in indistinguishability with coupling strength due 
to the Purcell effect, until it saturates for hg > 100 /zeV 
as the strong coupling regime is entered. The results from 
the ED and TCL clearly cannot be explained using this 
model. 

It is instructive to write the Hamiltonian, Eq. ([I|), in 
terms of polariton (dressed) states, which diagonalize the 
Jaynes-Cummings Hamiltonian, Eq. (l2|. This results 
in terms like <JiuJ2p ^^pi^p + ^p); which cause phonon- 
mediated transitions between the upper (u) and lower (1) 
polariton branches. These are separated by an energy 
Wu — wi = \/4:g'^ + A^, and hence this energy is expected 
to play an important role in the physical interpretation. 

Figure |2F a) shows examples of time-dependent dephas- 
ing rates calculated within the TCL, for two values of 
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FIG. 2. (a) Time-dependent dephasing rate from the TCL. 

(b) Effective phonon density at zero temperature. (c) 
Phonon population distribution function for hg = 30 /xeV, 
r = 0.5 ns-\ hK = 125 ^eV, and hA = 27.78 ^eV. (d) As 

(c) with hg = 200 ^eV. 



the coupling strength, g. Figure W[h) shows the corre- 
sponding effective phonon density, deflned as (T — 0) 
rfph(w) = '^J2p\^p\'^^i^ ~ ^p)- ^^^ both values of g, the 
dephasing rate attains large values within the first 3 ps, 
after which it settles to a smaller positive non-zero value. 
The initial temporal variations of the rate are directly 
related to a sampling of the entire effective phonon den- 
sity through virtual processes, which are allowed at short 
times due to the energy-time uncertainty relation. 

In the long-time limit, phonon- induced decoherence re- 
flects real phonon-mediated transitions, corresponding to 
the effective phonon density being sampled at speciflc en- 
ergies. This explains why the long-time value of the de- 
phasing rate is much larger for hg = 200 fieV compared 
to hg = 30 /xeV. The phonon density is thus sampled, 
respectively, at the energies given by the polariton tran- 
sitions, namely h{uja ~ ^i) = 2 x 30 /ieV = 0.06 meV, 
where the phonon density is small, and /i(aju — i^i) = 
2 X 200 /ieV = 0.4 meV, where it is much larger. The 
TCL only considers the long-time limit of the dephasing 
rate in Fig. ^a), explaining why the indistinguishability 
tends to unity for small QD-cavity coupling strengths. 

To verify the intuitive explanation provided by the po- 
lariton picture, we show in Figs, ^c) and (d) the phonon 
distribution function, {bj,{t)bp{t)) , calculated using the 
ED approach. For small QD-cavity coupling no speciflc 
phonon energy is singled out, consistent with the small 
phonon density at the corresponding energy of 0.06 meV, 
whereas for the larger QD-cavity coupling a signiflcant 
increase in phonon population occurs near 0.4 meV, as 
expected from the polariton interpretation. 

The nearly constant indistinguishability for the ED 
and IBM for small g arises due to the large difference in 



the involved timescales, the QD decay time, tqd = F 



off 



with 7 = 0, and the extent of the short-time non- 
Markovian regime, Tcorr [Fig. ^b)], where virtual pro- 
cesses dominate the decoherence [39] . The importance of 
the short-time regime is witnessed by the large dephas- 
ing rate in Fig. l2ja), which is especially important for 
the QD, and tqd needs to become comparable to Tcorr 
to affect the phonon dephasing. We note that for g -^ 
other phonon dephasing mechanisms can become impor- 
tant [nmS]. Further increasing g in Figs. Wic) and (d), 
real processes become increasingly important and con- 
tribute further to the decoherence, however, a stronger 
Purcell enhancement is also in effect, combating the in- 
fluence of decoherence by making the QD decay faster. 
Including only virtual phonon processes (IBM) , the Pur- 
cell effect can increase the indistinguishability, however, 
adding real processes (ED) the indistinguishability is seen 
to decrease monotonically, partly due the saturation of 
the Purcell enhancement, and hence tqd, in the strong 
coupling regime. Comparing the indistinguishabilitics for 
the QD and the cavity, only the ED predicts a significant 
difference between the two, which indicates that the dif- 
ference arises from short-time non-Markovian effects that 
are only retained in the ED. The smaller indistinguisha- 
bility found for the QD is a result of the direct interaction 
between the QD and the phonons, where the very strong 
short-time dephasing [Fig. [2]ja)] significantly decreases 
the indistinguishability of photons emitted from the QD. 
The photons in the cavity do not interact directly with 
the phonons, only indirectly through the QD-cavity in- 
teraction, and hence do not suffer to the same degree 
from the strong short-time dephasing as the QD does. 
Furthermore, the longer lifetime of the QD compared to 
the cavity, i.e. tqd ^^ 1/k, is also expected to have an 
influence as excitations residing in the QD simply have 
more time to interact with the phonons. We note that 
in the case of the phenomenological Lindblad theory, the 
same difference is not expected as here no short-time de- 
phasing is present, only the constant pure dephasing rate 
7- 

Spectral asymmetries. - Figures pTa) and (b) show the 
indistinguishability as a function of the QD-cavity detun- 
ing, A, which is an important experimentally controllable 
parameter. The Lindblad theory is unable to explain the 
variations with detuning that are predicted by the ED, 
both on a quantitive and qualitative level. The behavior 
of the Lindblad theory can again be understood using 
the analytical expression discussed above, since the Pur- 
cell enhancement decreases for increasing detuning. A 
common feature displayed by both the TCL and ED is a 
strong asymmetry with respect to the sign of the detun- 
ing. For large detuning, | A| » 5, the polariton dispersion 
becomes Wu/i ~ (A ± |A|)/2. Thus, to make real tran- 
sitions between the two polariton branches, the phonons 
need to provide an energy ±?i|A|, either through emission 
(— ) or absorption (-|-). At T = 0, only positive detun- 
ing will lead to phonon emisson and thus decoherence in 
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FIG. 3. (a), (b) Indistinguishability as a function of QD- 
cavity detuning for QD and cavity emission, (c), (d), and 
(e) Phonon population distribution function for different QD- 
cavity detuning. For A = the population has been scaled 
by a factor of 2. Parameters: KT — 1 fJ.eV, Hk = 100 /xeV, 
hg = 100 /ieV, and ^7 = 1.1 /leV. 



the long-time limit, hence a larger indistinguishability is 
expected for negative detuning. Indeed, both the TCL 
and ED display such an asymmetry. The effect of de- 
tuning on phonon emission is shown in Figs. P[c), (d), 
and (e), where, for fiA = 0.5 meV, a significant phonon 
population is observed at this phonon energy. 

While the detuning asymmetry of the cavity emission 
directly reflects the shape of the effective phonon den- 
sity in Fig. WCo)^ this is not the case for the QD emis- 
sion. The reason for this, perhaps surprising, difference 
between cavity and QD emission is the following: To gen- 
erate a cavity photon, the QD must decay by coupling 
to the cavity, i.e. subject to the Purcell effect. For large 
detuning (A ^ 5,k), the Purcell effect is only effective 
if assisted by phonon emission J9h11J. However, the QD 
can generate a photon without coupling to the cavity, 
namely through the background decay rate F. There- 
fore we expect the cavity to signiflcantly influence the 
QD emission only relatively close to resonance. For large 
detunings, we expect the QD result to converge towards 
the g — ;► result ^ , predicted in Fig. [ijc) . 

For negative detuning, the effect of virtual processes 
in the short-time regime is clearly seen in the ED re- 
sult, where despite the absence of phonon emission in 
the long-time limit, the indistinguishability is still signif- 
icantly below unity, especially for QD emission. This is 
not the case for the TCL which only describes the long- 
time limit. 

In conclusion, we have shown that non-Markovian 
phonon interactions strongly influence the coherence of 
single photons emitted from a cavity QED system. An 
exact diagonalization approach predicts an upper limit 
for the indistinguishablity, a feature not captured by the 
commonly used Lindblad theory. We provided physical 



insight into the non-Markovian dephasing processes us- 
ing a polariton picture. Finally, we predict an asymmetry 
in the indistinguishability with respect to the QD-cavity 
detuning. 
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